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Diseno de un generador eolico 



Analisis teorico del problema aerodinamico 

Para el presente trabajo se desarrollo una teoria sobre la mecanica de fluidos basada en un flujo 
incompresible, irrotacional y de viscosidad despreciable (no viscoso). Bajo estas hipotesis es valido el 
tratamiento del flujo de aire sobre las palas de la turbina como un flujo potencial. 

Definiciones y teoremas utilizados 

Se define la velocidad como una funcion (vectorial) dependiente del punto (x, y, z) y del tiempo. En 
nuestro caso suponemos el flujo como estacionario de forma que la derivada parcial de la velocidad respecto 
del tiempo es nula. 

V = V(x,y,z,t) 

dV 

— = — > flujo estacionario 
8t 

Considerando el flujo como incompresible se cumple que sobre una superficie de control el flujo del 
vector velocidad es nulo (entra la misma masa que sale del volumen de control). Aplicando el teorema de la 
divergencia, significa que la integral sobre ese volumen de control de la divergencia de la velocidad es nula. 
Como se puede elegir el volumen de control arbitrariamente significa que para que la integral sea nula debe 
ser nulo el integrando, de forma que la divergencia del vector velocidad es nula en todo punto. 

jv(x,y,z,t)-ds=0 

=> f div(v)-dv = 

JVol v > 

=> div^Vj=0 — > flujo incompresible 

Si el flujo considerado es irrotacional implica que la circulacion sobre una curva cerrada del vector 
velocidad debe ser nula. Aplicando el teorema del rotor significa que el rotor del vector velocidad es nulo 
para todo punto entonces la funcion de velocidad deriva de una funcion potencial de manera que: 

£v(x,y,z,y)-dl = 

Jrot(v)-ds = 

=> rot(v)=0 — >• flujo irrotacional 

rot(v)=VxV = 

^> V = VcE>(x,y,z)=cE> x i+cE> y ] + (i> z k 

div(v)=0 ^ div(vct>(x,y,z))=0 ^ V 2 cp(x,y,z) = 

Como el flujo es irrotacional e incompresible se desprende que la funcion potencial de la que se deriva 
la velocidad es armonica. 

Flujo bidimensional 

En nuestro caso se considera que el flujo en la coordenada z es nulo, de esta manera se puede tratar 
el flujo como bidimensional. 
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V = V(x,y) 

Escrita de esta forma se define la velocidad compleja como un numero complejo cuyas coordenadas 
son las componentes del vector velocidad. 

V = V(x,y) 

v = v 1 + i v 2 

Funcion Potencial complejo 

Como estamos considerando el flujo como irrotacional e incompresible podemos asegurar que existe 
una funcion potencial cuyo gradiente es la velocidad y ademas que la funcion es armonica. 

Entonces existe su armonica conjugada que representaran las lineas de flujo del potencial complejo. 

V = V(x,y) = VO = O x i+O y j 

V = O x +iO y 

V 2 ® = => $g Arm => B^P armonica conjugada de O 
=> 3 Q(z) = 0(x, y) + i ^(x, y) funcion potencial complejo 

Por la construccion de la funcion Q se puede demostrar que se cumplen las condiciones de Cauchy- 
Riemann (por ser ¥ conjugada armonica de <l>) entonces Q es una funcion holomorfa y se cumple: 

Q'(z) = © x +i^> x 

V = O x +iO y =O x -i^ x 

=> V = Q'( z ) 
Teorema de Blasius 

Al ser el flujo no viscoso se demuestra que cualquier obstaculo que se interponga en el fluido es una 
linea de corriente y la fuerza que ejerce el fluido sobre la superficie del objeto es siempre perpendicular a 
esta (la fuerza de rozamiento es nula). 

Sobre una linea de corriente se cumple la ecuacion de Bernoulli y si el flujo es irrotacional es valida en 
toda la masa de fluido. De estas hipotesis se desprende el teorema de Blasius que nos permite calcular la 
fuerza ejercida sobre un cuerpo interpuesto en un flujo irrotacional y no viscoso (por unidad de longitud) 
conocida su funcion potencial. 

V = Q'( z ) 

|dF| = p|dz| 
dF = i • p • dz 



F = i(tpdz 



Por la ecuacion de Bernoulli se tiene (sobre la linea de corriente y): 
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p + — pV =cte. 
2 ' ' 



=> F = i- 



cte 



iplvf 



V 



•dz = 



J 



-i -i 

= i • tt cte • dz — i • <p — plvl dz = -i- — pllvl dz 

Jy J?2 2 •'T 



dz = |dz|-e i(p 
F 



2 Jy e i2(|> 



dz-e 1. 



1 



-i- — p<fv 2 -Idzl-e^ 
2 h 

1 



F = i--p6V • dz -e ~ 1<p =i--p6V" -dz 

2 Jy I I 2 J Y 

Reemplazando la velocidad conjugada por la derivada del potencial se obtiene la ecuacion del teorema 
de Blasius. 



= l±pl(Cl'(z)Y-dz 



Transformation de Kutta-Joukovsky 

En este trabajo se utiliza una transformacion conforme que transforma una circunferencia desplazada 
del origen en un perfil aerodinamico llamada transformacion de Kutta-Joukovsky. Esta transformacion es 
conforme en todo punto salvo en el origen dado que la funcion no es holomorfa en ese punto. 

a 2 

W = ZH 
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Soluciones Particulares 

Torbellino Horario 



iy 
Q(z) = — lnz 

In 

Y i i Y Y i i 

Q(z) = z — (lnz + z'argz) = argz + z — lnz 

2.7t 2.7t 2.71 

Y 
Equipotenciales <p{x, y) = arg z = cte => arg z = cte 

In 

Y i i ii 

Lineas de Flujo \f/(x,y) = — lnz = cte => \z\ = cte 

In 




Cuerpo de Seccion Circular inmerso en un fluido con velocidad Vo(a = 0) 



w = z + - 



z = ae l( ? 

w = ae i( P + — = ae i( P + ae~ {( P = a.(e i(P +e~ iCP )= 2a 

ae l( P v ; 



cos cp 



V=Vo=Q'(w) 



Q(w) = Vo.w 



Q(z) = Vo\z + 
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Cuerpo de Seccion Circular (a = 0) + Torbellino de Aire 



( « 2 \ ; >i 

\z + — )+- ! —\nz 
\ z I 2n 



Q(z) = Vo 



Q.\z) = Vo\l- a - + r/ 



iS 



iS 



iS 



2n.z 



\ z 2 1 In.z 



.dz 



Vo 2 \l- 



+ 2Vo.l- 



a iy 



7 



z 2 )2n.z An 2 z 2 



dz 



Vo 2 1-4) dz + 2Vo.l^dz-l™^lz 

\ z 2 1 \2n.z \ 2n.z* 



An 2 z 2 



dz 



Vo 2 ll-\)dz = 



17 -dz = ^ 

2n 



2n.z 



2Vo.a 2 iy 

In.z" 



dz = 



4n 2 z 2 



dz = 



¥ = ^2VoS lr 



Serie de Laurent R z0 = 

iy 



Polo de l er Orden R, 



lim — z - 



2^- z^zo z 2n 

.2- 



Polo de 3 er Orden R = lim d Vo - a ^ i_ z 3 = o 



z ^ zo dz 2 ;r z ; 



Po/o de 2 ao Orden R = lim 



d x" 



'X 



dz = iS.Vo — = -iS.Voy 

2n.z 2n 



z ^ z0 dz 4n 2 z 2 
F = i.SVo.y.1 



: 2 =Q 




Cuerpo de Seccion Circular (a£0) inmerso en un fluido de velocidad Vo 
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z d) =ze ia 



Q(z (1) ) = -Voz (1) +4- -^lnz (1) 



Z W / ,271 



Q(z) = -voze ia + a ' ly -~- iot 



z e 



ia/ 2tc 



-lnze 1 



In z6 = In z + In 6 In 6 = cte desaparece al deri var 



Q( Z ) = -vo z e ia + a ' 1Y 



ze 



ia / 2tc 



lnz 




Perfil Aerodinamico 



Se define entonces una circunferencia desplazada y aplicamos la transformacion. 

V 
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R. 



•Va 2 +m 2 



R = Va 2 +m 2 +d 

^ = z + — 



Puede verse que la circunferencia de radio Re se transforma en el "esqueleto" del perfil aerodinamico 
y se puede demostrar que es un arco de circunferencia. La circunferencia exterior de radio R se transforma 
en el perfil aerodinamico. 

Segun el teorema de Blasius la fuerza ejercida sobre el perfil sera: 

^ r 
Para poder calcular esta integral es necesario un cambio de variable: 

>* Izl 




z = z (1) +d.cosP + i(n + d.senp) 
z (1) =z-d.cosp-i(n + d.senp) 

^ = z + — 
z 

De donde es conocido el potencial complejo para la variable z 1 asi como la velocidad expresada como 

derivadas de funciones compuestas. 



Q(z w ) = -Voz w e ux + 



\ z (1) e ia ' 2n 



lnz 



(i) 



dQ dQ dz (1) dz 



V: 



dQ „ Li„ R 2 



dz (1) 

dz (1) 
dz 



d^ dz (1) dz d^ 

■=-vo(e ia 



iy 



= 1 



z (D 2 e ia/ 27i.z (1) 
z (1) =z-d.cosp-i(n + d.senp) 
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Aqui es necesario destacar que la velocidad debe ser finita en todo punto, analizando los puntos 
singulares de la derivada de z respecto de % se desprende que tiende a infinito cuando z tiende a (-a). Para 
salvar este problema y que la velocidad sea finita en ese punto la derivada del potencial complejo respecto 
de la variable z 1 debe anularse en ese punto. De esta manera se puede calcular facilmente el valor de la 
circulacion que genera el perfil (y). 



Q'(-R.e 1 P) = -Vo 



/ 



»ia 



R' 



\ 



V 



vo e 



iia _ p -ia p -2ip 



^2p2iPaia 

iy 



ly 



J 



271.R.6 



iP 







dz 



d^ 



->00 



r.vo e^e 



27t.R.e 

iapiP p-ia P -ip^ = ' l 1 

2n 



ip 



y = 2tt:.R.Vo 



' e i(a + p)_ e -i(a + p)' 

y = 47t.R.Vo.sen(a + p) 

Entonces es posible calcular la integral de Blasius para averiguar la fuerza ejercida por el perfil. 



r 
F = ^|(Q'(z)) 2 ^'(z)-dz 



F = ^(n'(z)r. 



( ~2\ 



.dz 



J 



F = ^. i (Q- ( Z )) 2 dz = i.8.Vo.y.l 
2 r 



V Z J 



no aporta Re siduos 



F = p.Vo.y.l 

F = p.Vo.l.47r.R.Vo.sen(a + p) 

F = 4.7i.p.Vo 2 .l(Va 2 +n 2 + d)sen(a + p) 



F = 4.7t.p.Vo 2 .a.l. 



1+ 



'nY d 



w 



sen(ot + P) 



F = 4.7i.p.Vo 2 .a.l. Jl+(tgp) 2 +- lsen(a + p) 

V a) 



Consideraciones sobre potencia del viento 

Teorema de Betz 

Para calcular las dimensiones que debe tener la turbina es necesario saber que cantidad de viento 
debe pasar por ella y cuanto de este es aprovechable como energfa mecanica. Para ello se considera la 
potencia del viento que fluye por un cilindro de radio 'D' siendo su area: 



Pagina 8 de 33 



5H FACULTAD 
5* DE INGENIERIA 

■^h Unive^sifJad de Buenos Aires 



Joaquin Ignacio Aramendia 

Padron: 87651 

e-mail: samsagax@hotmail.com 



A = n- 



^ 



La energfa de ese tubo de viento es: 
E=-m-V 2 



Ysu potencia: 

p= E = lm 

t 2 t 



V 



■V 



•V' =A- 



2 



Considerando la interposicion de la turbina en el tunel de viento se pueden considerar cuatro 
secciones, una muy alejada de la turbina donde la velocidad del viento se considera igual a la que existe si la 
turbina no estuviese (1), una muy alejada de la turbina (4), pero despues de esta donde la velocidad sera 
modificada pero la presion sera la misma que en la primer seccion (P a tm)- Una seccion justo antes de la 
turbina (2) y una justo despues (3) donde no se puede inferir que pasa en el medio pues se le extrae una 
cierta potencia al tubo de viento. 




De esta manera se puede aplicar la ecuacion de Bernoulli entre las secciones 1 y 2 y entre las 
secciones 3 y 4 (pero no entre las secciones 2 y 3): 



2 

1 



2 

1 



P 3 +^pV 3 2 =P 4 +-pV 4 2 

Siendo el flujo masico en todas las secciones el mismo y que la presion estatica de la seccion 1 es la 
misma que la seccion 4 

P 1 =P 4 =P 

A lP V x = A 2 pV 2 =A 3 pV 3 = A 4 pV 4 



A 2 = A 3 
1 



V, s V, = V 



1 



P + -pV, --pV / 



1 



1 



p + -pV — pv 



p -P 



ipK 1 



V/) 
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P VA(v 1 -V 4 )=F = ipA(v i 2 -V 4 2 ) => V = &±^ 



De aqui se puede calcular la potencia util 

Pot u =F-V = ^pA(v i 2 -V / 2 

Que comparando con el resultado anterior de potencia se puede definir un rendimiento: 



IpA(v=-V/)&±^ 



a 



n 



Pot u 
Pot 



ipAV^ 



1- 




2 


r v 4 ^ 



[l-a 2 ](l+a) 



Para averiguar el rendimiento maximo se deriva la expresion respecto de la variable 'a' e igualando a 



cero: 



dr| _ -2a(l+a) + l-a 2 



da 



a = 



Y con este valor se saca el rendimiento maximo, llamado rendimiento de Betz, utilizado para 
dimensionar la turbina. 



11 



1- 



fi\ 



V-V 



Ml 



16 

27 



= 0.592 
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Analisis teorico del problema estructural 

En el presente trabajo se calculara estructuralmente la turbina basandose en un analisis estatico y un 
analisis dinamico para calcular las frecuencias naturales de oscilacion de las partes mas importantes del 
generador. Para todo el desarrollo se toma como hipotesis que las deformaciones son elasticas, donde es 
valida la ley de Hooke. 

Problema estatico 

Viga sometida a esfuerzo axil 







Al 










F 



























Dada una configuracion de una viga empotrada en un extremo y cargada en el otro con una carga 
segun su eje se define la tension en cada seccion como la carga dividida el area de la seccion. 

F 

a = — 

A 

Esta carga dara lugar a una deformacion que se puede escribir en terminos especificos dividiendo la 

deformacion total por la longitud de la viga, a este cociente se lo denomina alargamiento especifico o 

deformacion especifica: 

Al 
s = — 
1 

Existe una proporcionalidad entre la deformacion especifica y la tension aplicada dada por una 

constante del material (su modulo de elasticidad, o modulo de Young) que da lugar a la ley de Hooke: 

a = E -8 

Viga sometida a flexion 




Dada una viga con una carga en funcion de su longitud se define el esfuerzo de corte como la integral 
de la carga distribuida: 
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Q = }q(x)dx 



Y el momenta flexor como la integral del esfuerzo de corte: 

X 

M = j*Q(x)dx 



Bajo la hipotesis de que en la flexion la seccion se mantiene plana y gira respecto de un eje llamado 
neutro (donde no existe tension) se asimila el perfil de tensiones sobre una seccion como un perfil lineal: 

y M 

Entonces se tiene para un diferencial de area de seccion un diferencial de fuerza, integrando ese 
diferencial de fuerza multiplicado por la distancia al eje neutro, se obtiene el momenta flexor en esa seccion, 
lo que permite relacionar el momenta en la seccion con la tension maxima aplicada: 

dF = a • dA 



dM = dF • y = a • dA • y = -^ y'dA 

y M 



y M 



J A yMA 



cu 



Donde J es el momenta de inercia de la seccion respecto de la linea del eje neutro (perpendicular al 
piano de aplicacion del momenta flexor). De esta forma se obtiene una expresion para verificar la seccion a 
flexion conocidos su momenta flexor y su momenta de inercia, lo que nos permite dimensionar la seccion. 

M 



CT, 



Curva elastica de una viga en flexion 

Conocidas las deformaciones especificas maximas en las distintas secciones de la viga es posible 
determinar el angulo de curvatura de la viga. Considerando que las secciones se mantienen planas luego de 
deformadas, el angulo de curvatura sera el determinado por dos secciones de la viga. 
a, 



Adx 



'M 

E 



a M dx M y r 



'M 

E 



E-Jy t 



dx = 



M 
E-J 



dx 



dcp = 

y M 

cp = ar tan(y') 

dcp_ y" 

dx 1+y' 2 

Considerando las deformaciones muy pequehas se puede despreciar la derivada primera de y al estar 
elevada al cuadrado frente a 1, esto permite aproximar el angulo a la tangente entonces es posible escribir: 

<p = y' 

Y por las relaciones anteriores se tiene: 
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y. 1 1 we ■ sic : ; 1 1 1 de Buenos Aires 


d 2 y 


M 


dx 2 


E-J 


d 3 y 


Q 


dx 3 


E-J 


d 4 y 


_q(x) 



dx 4 E-J 

Ecuaciones que nos permiten, integrando, conocer la curva elastica de la viga conocida la carga 
distribuida sobre esta. 

Problema dinamico (vibraciones en solidos) 

El problema dinamico consiste en determinar las frecuencias de oscilacion naturales de cada elemento 
estructural basico del generador (palas y columna) a fin de evitar que el funcionamiento normal del mismo 
lo Neve a una falla. 

Es preciso determinar las frecuencias de oscilacion axiales en la columna para que no entre en 
resonancia debido a la velocidad de giro de le helice y algun desbalance estatico en el centra de mansa de la 
misma. Asi mismo es necesario determinar las vibraciones de flexion en la columna con los mismos fines. En 
las palas es preciso hacer un analisis de las frecuencias naturales de vibracion en flexion, dado que cada vez 
que una pala se enfrenta con la columna del generador produce una reduccion en la sustentacion que Neva 
a que el momento fexor caiga subitamente. 

Vibraciones axiles no forzadas (frecuencias naturales) 

Cuando la carga sobre el extremo de la viga empotrada varia en el tiempo la deformacion sera funcion 
no solo de la posicion sino tambien del tiempo por la ley de Hooke. Llamando Y a la funcion de 
desplazamiento para una posicion V a un tiempo Y, entonces la variacion de Y esta determinada por: 

y = y(x,t) 

Ay = y(x + Ax,t)-y(x,t) 

Siendo la deformacion especifica Y el cociente de la deformacion y la longitud inicial: 

_ Ay _ Ay = y(x + Ax,t)-y(x,t) 



Ax 




Ax 






Ax->0 










dy 
s= — 

dx 










Entonces la tension aplicada 

a = E -8 = E^^ 

dx 


en una seccion 


V sera 



Analizando la fuerza sobre el mismo elemento de barra 'Ax': 

F 1= a(x,t)-A 

F 2 = a(x + Ax, t) • A 

Donde A es el area de la seccion de la barra. Por la segunda ley de Newton se tiene: 
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8 y 

F x = m • — J - 

dt 2 



d 2 v 
A ■ [a(x + Ax,t)- a(x,t)] = A • Ax • p • — y - 

di 



AE- 

dy 



^(x + Ax,t)-|(x,t) 

<3x ox 



= A • Ax • p 



dt 2 



fy , 



. (x + Ax,t)--^(x,t) 

<3x ax 



Ax 
Ax^O 

d 2 y = p d 2 y 
5x 2 E ' dt 2 
^y _ E 5 2 y 



p d 2 y 

E St 2 



p 8x' 



c 2 = 



Entonces la funcion de desplazamiento corresponde a una ecuacion de onda unidimensional, en este 
caso longitudinal. La resolucion de esta ecuacion es posible por separacion de variables de la siguiente 
forma, planteando una solucion como el producto de dos funciones de una sola variable: 

y = x(x)-T(t) 

T"X = c 2 X"T 

Sin perdida de generalidad se aplicaran las siguientes condiciones de contorno: 

y(0,t)=0 

a(l,t) = 0=>y(l,t)=0 

y,M) = o 

y(x,0) = ^x 

Lo que da lugar a un problema de Sturm-Liuville en la variable V: 
y(0,t)=X (0) T (t) =0=>X (o) = 

y(U) = x (1) T (t) = o^x (1) = o 



y,M) = x (x) T ( ' 0) = o^T/ 0) = o 



y(x,0) = |x 

Separando las variables: 
T"X = c 2 X"T 

:r _^L_r _ 

X ~c 2 T ~ 

jx"+xx = o 

[T" + Xc 2 T = 

Cuya solucion en X es la combinacion lineal de exponenciales que puede escribirse en forma de senos 
y cosenos: 
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X„ = 



'n^ 2 



vly 



Siendo % X, n ' los autovalores del problema. La funcion T' tendra una forma similar, de lo que se puede 
escribir a Y como una serie de Fourier. 



T(x) = C sin(VXct)+ D cos(V^ct) 
T'(o) = C = 



y( x >On =B n sin 

CO 

y(x,t) = ^B n sin 

n=0 

CO 

y(x,0) = XB n sin 



— x 

vl J 

r 



•cos 



nn 



-ct 



n=0 
1 



rm 
T 
rm 
T 



V 



cos 



rm 



-ct 



Al 



I 



Al 



V J- J 
f 



x -sin 



=> B„ =-5- 



V 



rm 

T 



■dx 



sin 



'rm 
— x 
1 



V 



2 rAl 
— x-sin 



'rm > 



l J o 1 



•dx 



Sin embargo no es necesario determinar los valores de los coeficientes pues las frecuencias naturales 
de oscilacion se pueden sacar del termino que multiplica a la variable Y 

rm E 



co r 



rm 
— c = — 
1 1 



Vibraciones axiales de una columna con una masa concentrada (frecuencias naturales) 

En nuestro caso deben considerarse las vibraciones naturales de la columna con la masa de toda la 
maquinaria del rotor supuesta concentrada en el extremo. Es preciso que las frecuencias naturales no se 
encuentren en el rango de giro de la helice de la turbina. 

Se supone toda la masa de la helice y de generador concentrada en el extremo de la columna como 
una masa 'm'. En este caso se tiene las siguientes condiciones: 

F = Q-A = E-A-y x (l,t) = -m-y tt (l,t) 

y(0,t)=0 

y,(x,o)=o 

y( x >0= x ( x ) T (t) 



ry>ff 


X = c 2 - 


X' 


'■T 


X" 
X 


1 T" 


= 


-X 


f x 


' + XX = 







1r 


+ Xc 2 T 


— 
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x(x) = Cj sin(vXxj+ C 2 cos(yXxJ 

x(o) = c 2 =0 

T(t) = C 3 sin(V^ct)+ C 4 cos(V^ct) 
T'(0) = C 3 =0 

X(x) = Cj sin(yXxJ 

T(t) = C 4 cos(VXct) 
E-A.y x (l,t) = -m.y tt (l,t) 

E-A-X^.T^-m-X^T^ 

E • A-X/a •C 4 cos(yXctJ=m-X( 1 ) -A,-c 2 -C^os^ActJ 
E • A • yfk ■ C i cos(yAAj = m • C 1 sinU/Xlj- X ■ c 2 
E • A • cos(Vxi) = m • sin(Vxij- vT • c 2 



EA Ap 



X tan 



(V5ii) 



m-c m 

Si llamamos 'B' a la masa de la barra y V al argumento de la tangente: 

l-Ap = B 

Jx-\- 

B 



m 



a 
a • tan(a) 



De donde se pueden resolver numericamente los valores de 'a' para obtener los valores de las 
frecuencias naturales de oscilacion del sistema. 

Vibraciones axiales forzadas 

Como debido a la accion de una masa excentrica girando a una cierta velocidad angular se genera una 
fuerza pulsante de la misma pulsacion (la frecuencia de giro) se da el caso de vibraciones axiales forzadas. 
Es necesario que las frecuencias tales no sean coincidentes con las frecuencias naturales para evitar que la 
estructura entre en resonancia. 

La ecuacion diferencial a resolver es la misma pero se cambian las condiciones de borde, donde esta 
vez la tension sera una cierta funcion del tiempo. 

y(o,t) = o 

F • sin(cot) _ / n, F • sin(cot) 



A 



y x (u): 



A-E 



a(l,t): 

y t (x,o) = o 
y(x,o) = fx 

Nuevamente podemos expresar la funcion Y como el producto de dos funciones de una variable para 
separarlas. 

y = x(x)-T(t) 

T"X = c 2 X"T 



= -X 




T' + Xc z T = 
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Cuyas soluciones son conocidas pues es el mismo sistema con otras condiciones de borde. 
y(0,t)=X (o) T (t) = 0=>X (o) = 

x(x)= C 1 sin^vA-xj+C 2 cos^A-x] 

x(o) = c 2 =0 

T(t) = C 3 sin(VIct)+ C 4 cos(V^ct) 

y(x, t) = sin^vAxJC 3 sin^vA,ctj+ C 4 cos^/ActJ 

y x (x, t) = VA- cos(yA,x|C 3 sin(yA.ct]+ C 4 cos(yA,ctj] 

y x (l,t)= V^ cos(V^l|c 3 sin(V^ct)+ C 4 cos(V^ct)]= ^Hgj => c 4 = 

A • E 

sin(oot) = sin^VA-ctj => vAc = co 

= C,VA cos(vAl) =^> C, = = 1 _\ 

A-E V ' A-E-VIcos(Vai) 

^ y ( x ' ^ = a ^ rv — nrn sin ^ x )' sin (^ ct ) 

y (l,t)= = 1 = \ cos(-\/Ax)-sin(->Act)= -t = \ cos(vAx)-sin(VAct) 

A-E-VAcos(Va1J V ; V ; A-E-cos(Va1J V ; V ; 

De esta forma podemos determinar la tension dinamica del punto mas solicitado que sera el de la 
base de la columna. 

y x (x, t) = f-^r^cosUXx )• sin(oot) 

A-E-cos(VAlJ 

a = E • y x (0, t) = p-^ • sin(cflt) 

A • cos^v XI) 



max 



COS 



(S) 



La fuerza F debe considerar la masa total de la helice mas el generador y toda la maquinaria. De 
forma que una vez calculadas las vibraciones naturales de la columna con la masa total en el extremo libre, 
es necesario que la frecuencia de giro del rotor este suficientemente alejada (mas de un 10%) de la 
frecuencia de resonancia del sistema. 

Vibraciones de flexion (frecuencias naturales) 

Se calcularan las vibraciones de flexion tanto en las palas como en la columna para elegir una 
velocidad de rotacion suficientemente alejada de las frecuencias naturales de oscilacion. En ambos casos se 
tomara el momento de inercia de la seccion como constante, en los casos en que no lo sea (como sucede en 
las palas) se tomara el momento de inercia de la seccion media como constante. 

Considerese una carga distribuida sobre una viga empotrada en un extremo donde la ecuacion 
estatica es la ya deducida anteriormente. Con un razonamiento similar al de las vibraciones axiales, 
desarrollando la segunda ley de Newton se tiene: 

E,.fX qM -A.,ff 

Donde en este caso la densidad es por unidad de longitud. La carga aqui puesta como una funcion de 
la posicion y del tiempo es en realidad una carga que fuerza vibraciones. Como nuestro caso de estudio es 
de vibraciones no forzadas esta carga sera nula de forma que nos queda una ecuacion mas simple: 
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5V 
dt 2 



Ap 



y xxxx y tt 

En el caso de una viga empotrada las condiciones de borde son las siguientes: 

y(o,t) = o 
y x (o,t) = o 
y xx (u) = o 
y xxx (lt)=o 

y(x,0) = f(x) 

y t M)=o 

Las primeras cuatro condiciones corresponden al estado de cargas y de vinculo. Las ultimas dos 
condiciones se explicitan pero no influyen en el calculo de el calculo de las frecuencias naturales de 
oscilacion. 

y = X(x)-T(t) 



a -T-X" 

, X"" 



X 



CO 



= -XT" 

— = X = co 2 
T 



X = 



X" 

a z 

T" + co 2 T = 

Analizando las condiciones de borde del problema: 



X 



(o) 







y(0,t) = X (o) T (t) =0 

y,(o,t)=x; 0) T (t) = o 
y„(i,t)=x; ) T (t) = o 

y„(l,t)=X5)T (t) = 

Se propone una solucion exponencial que desarrollada queda: 
TV t ) = Acos(cot) + Bsin(cot) 

T/ t ) =-coAsin(cot)+coBcos(cot) 



X' (0 ) = o 

x w = o 



V, 



(o) 
r (0 



coB = 
B sin (cot) 



B = 
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X = e rx 
X"" = rV x 

UJ 



r 4 e rx 



f-1 



= 



r = *i 



'<d^ 



vay 



X = Ccos 



X (o) =C + E = 



A 



J 



+ Dsin 



f \ — A 
CO 
X 

a 



v 





'aA 


i 


vaj 


J 


^ 


I 




vaj 


<• > 






' (i) , 




J 


— 


•l 


vU; 


J 


'co^ 


. V 


vaj 


+ 


Ecosh 



X' 



Esin 



f i — A 

CO 

•x 
a 



C = -E 
+ Dcos 



V 



r i — A 

CO 

•x 
a 



f i — A 

CO 

x 
a 



v 



+ Fsinh 



f i — "\ 

CO 

x 



\ 



+ E sinh 



f \ — A 
CO 
•X 

a 



+ Fcosh 



/ 


(CO 


\ 




- 


• X 





la 


y_ 



X' (o) =^[D + F] = => D = -F 



X" 



x '('i) 



X'" = 



CO 



cos 



/i — \ 

CO 
•X 

a 



+ cosh 



v 



CO 



'coV2 



u; 



cos 



J 

f i — A 
CO 

v .a j 

f 
sinh 



f i — A 

CO 
•X 

a 



+ cosh 



v ,a j 



+ F 



+ F 



sin 



f I — A 
CO 

•x 
a 



sin 



v 

^ i — \ 

CO 

v ,a j 



+ sinh 



f i — A 

CO 

x 
a 



+ sinh 



sin 



+ F 



cos 



f i — A 
CO 

X 

v .a j 



V 

f I — -• 

CO 

v .a j 

r 
+ cosh 

v 



Y'" 



'coV 2 



w 



sinh 



■1 



sin 



v .a , 



■1 



v ,a j 



+ F 



cos 



^ I — \ 

CO 



+ cosh 



^ i — A 

CO 



= 



De donde se puede armar el siguiente sistema para % E' y T', cuyo determinante debe ser cero: 
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cos 



sinh 




+ cosh 



-sin 




+ sinh 



+ cosh 





[cos(u) + cosh(u)] - |sin 2 (u)+sinh 2 (u)J = 



cos(u) + cosh(u) sin(u) + sinh(u) 
sinh(u)-sin(u) cos(u) + cosh(u) 
cosh 2 (u)+ 2cosh(u)cos(u)+ cos 2 (u)- sinh 2 (u)+ sin 2 (u) = 
cosh(u)cos(u) = -1 

De donde, despejando para V se pueden sacar las frecuencias naturales de oscilacion V. La 
velocidad de rotacion debe estar suficientemente alejada de esa frecuencia (preferentemente por debajo de 
la primera armonica) para evitar la resonancia y posterior rotura de las palas o la columna. 

Es posible determinar la solucion completa para y despejando 'E' en funcion de T' y reemplazandola 
en la ecuacion de y, lo que daria una forma de serie de Fourier. Conocida la funcion inicial de la posicion de 
la viga, es posible determinar los coeficientes de la serie, sin embargo no es el proposito de este estudio la 
funcion en si, sino solamente las frecuencias de resonancia de la viga. 
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Analisis practico del problema aerodinamico 

En esta seccion se aborda el tema de la deduccion y fijacion de parametros de diseno del perfil de 
cada pala de la turbina eolica. El objetivo de este trabajo es dimensionar una turbina eolica que genere 3 
Mw de potencia electrica considerando el rendimiento de Betz y el rendimiento de la maquinaria involucrada. 
En una seccion posterior se analizara el problema estructural considerando no solo el analisis estatico sino 
tambien el problema de vibraciones naturales de la estructura. 

Para determinar las condiciones de funcionamiento y determinar las dimensiones principales de la 
turbina se elegira una ubicacion. En este caso la ubicacion se eligio en la zona oeste de la ciudad de 
Comodoro Rivadavia, donde los vientos son los mas fuertes de la Argentina, con promedios anuales de 11,2 
m/s. 

Potencia mecanica requerida por el rotor y radio de barrido de palas 

Considerando la potencia electrica nominal (Pe), el rendimiento de la maquinaria del generador y el 
rendimiento de Betz es posible calcular la potencia mecanica que debe generar el rotor asi como el area de 
barrido de las paletas y por tanto, su longitud. Del teorema de Betz tenemos: 

TW =^ = 0,592 

v 

P^p-A-V, 3 



A 



2 • P„ 2 • P., 



P-V/ p-V^-r)^ 




Considerando el rendimiento de la maquinaria cercano al 0.85 
P„ =3Mw = 3-10 6 w 




2P„ 

= 47.2257m 



"•p-V 1 3 -r )Betz -ri r 



I maq 

De este calculo adoptamos un radio de pala de 48m por simplicidad en los calculos, tambien es 
necesario adoptar un radio minimo correspondiente a una porcion inactiva de la pala, para este caso se 
adopto 5m para evitar todo tipo de interferencia entre las palas. De todas formas las fluctuaciones del viento 
hacen que la potencia entregada no sea constante. La forma de controlar este aspecto es con el angulo de 
ataque de las palas respecto al viento variable en funcion de la velocidad del mismo. 

Se busco en turbinas parecidas las caracterfsticas usuales de las palas para definir las dimensiones. 
Adoptado un generador sincronico, de manera que la velocidad de giro del rotor se mantiene constante, es 
necesario reducir la velocidad de rotacion utilizando una caja reductora de engranajes, un freno mecanico 
que se activa en condiciones adversas (si es necesario parar la turbina por vientos demasiado fuertes, 
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mantenimiento etc.) y un regulador del angulo de ataque de las palas para regular la potencia. La velocidad 
del rotor se fija en 15 rpm, lo que hace necesaria una caja de engranajes que multiplique esa velocidad por 
100 para un generador de 4 polos magneticos. 

Para el ancho de palas se propondra un perfil lineal con el radio (el parametro 'a' sera lineal, lo que 
aproximadamente controla el ancho del perfil). La relacion \j/a' se mantendra constante a lo largo de toda la 
pala y se toma entre 0.1 y 0.2 para no hacer que el espesor del perfil sea excesivo (lo que aumenta la 
fuerza de 'drag' que resulta contraproducente) y tampoco muy pequeho que comprometa la rigidez 
estructural de las palas. De la investigacion realizada se adopto una pala de ancho 3m en el radio minimo y 
de lm en el radio maximo (a max = 0.75 y a min = 0.25 respectivamente a los radios R min = 5m y R max = 48m). 

Curvatura de las palas 

Para analizar la curvatura de las palas es necesario tener en cuenta que la velocidad aparente del 
viento sobre la pala es la composicion vectorial de la velocidad de rotacion de la pala y del viento incidente 
sobre esta, de manera que la velocidad aparente sera variable sobre el perfil y es necesario que la pala 
tenga un angulo de posicion respecto del eje del rotor variable con el radio. Para eso definimos el angulo V 
como el angulo formado por la velocidad aparente y el eje del rotor. Se adoptara el angulo 'p' (curvatura del 
esqueleto del perfil) y 'a' (angulo de ataque del perfil) como constantes para el diseno. 



I-JH J »l I I 




V rot =co-R 
V„ =11.2 m / 



cte. 



V a=A /(co-R) 2 +V v 2 

De manera que el angulo de rotacion de cada perfil (definido como el angulo que formara el eje del 
perfil con el eje del rotor) en funcion del radio estara dado por 'X' (funcion del radio) y V (tornado constante 
en cada seccion). 



^(r) 



X( R ) +a = atan 



CD-R 



V. 



+ a 



v J 



Esto hace que la fuerza sea variable tambien con el radio a medida que se aleja del eje del rotor. 

Momento generado por las palas 

La fuerza util en la generacion de potencia es la componente perpendicular al eje del rotor, como la 
fuerza resultante es perpendicular a la direccion del viento aparente se tiene: 
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X 



(R) 



a tan 



co-R 
V 



F u =F-cos(x (R) )=F- 
Fper =F-sin(x. (R) )=F- 



co-R 



V., 



Notese que la direccion de la fuerza esta determinada solamente por el angulo que forman la 
velocidad aparente del viento y el eje de la turbina. De la parte teorica tenemos una expresion para la fuerza 
por unidad de longitud de perfil, de la que integrando en 'R' la fuerza util es posible determinar el momento 
util que genera cada pala, multiplicado por la cantidad de palas (en nuestro caso 3) se obtiene el momento 
total que genera el rotor: 

r A , ^ 



F = 4-7t-p-a- - + -y/l+tan 2 p • V a 2 -sin(a + p) 



Va 



K = 0.15 



a =b-R + c 



R — R 

max min 

C = a max - t>-R min 



-\ 



F u = F^- = 4 • 7i • p • (b • R + c)- - + ^/l+tan'p 1 • V a • V v • sin(a + p) 

v a U J 



R m » R max /J \ 

M u =3- JF U -R-dR = 3- j 4- rc-p-(b -R + c)- - + 7l+tan 2 p -V a • V v -sin(a + p)-R -dR 

R . R . V 3 ) 



M u =12-TC-p-a- - + A /l+tan 2 p • V v -sin(a + p)- [l (Rmax) -I (Rmin) ] 



Va 



Donde: 



I (R) =J(b-R + c\/(co-R) 2 +V v 2 -R-dR 



fy ^ 



= |o)-b-R 2 - (R) 2 + -*- dR + jco-c-R- (r) 2 + 



V COy 



V, 



CO 



dR 



8 = ^ 

CO 



I 



cob 



(R) 



R(R 2 + 8 2 )^ + ^(R 2 + 8 2 )^-^ln[R + (R 2 + 8 2 )^ 



+ ^(r 2 +S 2 )^ 
3 v ; 



Con esta expresion se puede calcular la potencia mecanica generada por el rotor, que teniendo en 
cuenta el rendimiento de la maquinaria (estimada en 0.85) daria la potencia electrica generada. Se utilizo 
una funcion programada en Matlab que se adjunta en el Apendice. 

Para los parametros fijados, la potencia electrica sera: 
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a = 5° = cte. 
p = 5°=cte. 



R max =48m 



R. 



5m 



a max = 0.75m 



a max = 0.25m 
V =11.2 m / 



oo = 15rpm 



K = 0.15 



TCl 

2 s 



M u = 2.3662 -10 b Nm = 2.3662 MNm 
p u =M u -oo = 3.7167 Mw 
P e =P u -ri maq =3.1592Mw 



Valor de Pe superior al buscado. Notese que el efecto del angulo 'a' es predominate sobre la potencia 
electrica. En efecto si se cambia el angulo a 4° la potencia electrica cae a 2.8461 Mw. De manera que 
equipar a la turbina con un sistema que pueda girar las palas para aumentar o disminuir el angulo de ataque 
en funcion del viento resulta muy ventajoso; ademas es una de las formas mas efectivas de controlar la 
potencia electrica entregada por el generador, actuando sobre la fuerza que generan las palas sin variar la 
velocidad del rotor. 

Modelo constructive de una pala 

Una vez fijados los valores de los parametros del perfil, es necesario llevar a cabo el modelo de una 
pala. Para ello se efectuo un programa que dibuje el perfil (adjunto en el Apendice), posteriormente se le 
agrego un modulo para que dibuje perfiles en funcion de la distancia al rotor y le de a cada seccion de la 
pala la rotacion de acuerdo al viento aparente (manteniendo el angulo 'a') constante en todas las secciones. 
Luego de obtener varios perfiles se pasaron los resultados a un programa de CAD para obtener un modelo 
en 3d de una pala. 

Usando la transformacion de Kutta-Joukovsky con los parametros fijados se llego a cada perfil, aqui se 
muestran el perfil generado para el radio minimo y el radio maximo sin su rotacion correspondiente por 
claridad. 



Perfil 
R=5m 
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Perfil 
R=48m 
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Con el mismo programa se genero una vista preeliminar de las distintas secciones (separadas cada 
lm) con su respectiva rotacion. En el Apendice se adjunta ademas el modelo 3d del esqueleto de la pala 
propuesto para resolver el problema estructural. 



Pala 3D 




x [mm] 
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Analisis practico del problema estructural 

El problema estructural se reduce a un analisis estatico de las palas y la columna y el analisis dinamico 
en los distintos modos de solicitacion de los mismos elementos. 

Para reducir el peso de la pala se opto por no fabricarla maciza sino tomando a intervalos regulares 
chapas con el perfil disenado uniendolos con dos almas longitudinales. (Se adjunta un modelo en CAD y los 
datos extrafdos del mismo en las distintas secciones en el Apendice) 




Diseno de pala 
Analisis estatico de la pala 

Con la expresion que se tiene de la fuerza total, es posible calcular el momento flector integrando dos 
veces la fuerza distribuida en funcion del radio. 

F = 4 • 7i • p • (bR + c)- (k + Vl+tan 2 p)- V a 2 • sin(a + p) 

Q(R) = JFdR = 4 • 7i • p • (k + ^/l+tan'p)- sin(a + p) • j" (bR + c)((coR) 2 + V v 2 )dR 

q(r) = 4 • n • p • (k + Vl+tan 2 pj- sin(a + p)- [l 1 (R)+ cj 
I 1 (R)= j(bR 3 co 2 +bV v 2 R + cR 2 ro 2 +cV v 2 )-dR = 

I 1 (R) = ^R 4 + ^R 2 + ^R 3 + cV v 2 R 
v ' 4 2 3 

M( 



(r) = JQ(R)dR = 4 • 7i • p • (k + Vl+tan 2 pj- sin(a + p) • J [^(r) + C x ]dR 



M(R)=4-7i-p-(K + Vl+tan 2 p)-sin(a + p)-[l 2 (R) + C 2 ] 



I 2 (R)=f[l 1 (R) + C 1 ]dR = ^R 5 + ^R 3 + ^R 4 + ^R 2 +C 1 R 

2V / JL1V / U 2Q g 12 2 

Con las condiciones de cargas en el extremo nulas: 

Q(R max ) = o => [i 1 (R max )+c 1 ] = o ^ Cr-^R.) 

M(R max ) = => [l 2 (R max )+C 2 ] = ^ C 2 =-I 2 (R max ) 

La seccion correspondiente al radio minimo (la mas solicitada) debe ser verificada a flexion, siendo el 
momento de inercia el sacado del programa de CAD y el modulo de Young correspondiente al acero. Se 
puede entonces calcular la tension maxima en esa seccion. Se considerara el peor de los casos que 
corresponde al minimo momento de inercia de la seccion, dado que el momento a lo sumo se encontrara 
actuando en ese piano. 
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R 255,973 mm 




Vista de section R=5m 



M (R ml J 

_ mm/ 

max T J r 



j = I x = 0.00910887 m 
y max = 0.225973 m 



M (RmJ = 2.8397 -10 b Nm = 2.8397 MNm 
a max = 70.44733 MPa 



Curva elastica de la pala 



Considerando el Momenta de inercia de la seccion a lo largo de la pala como constante e igual a un 
momenta de inercia medio (el valor es tornado de la seccion media de la pala, cuyas caracterfsticas se 
adjuntan en el Apendice) es posible determinar la cuerva elastica de la pala. Nuevamente tomamos el peor 
de los casos en el que el momenta flector actua sobre el piano correspondiente al momenta de inercia 
minimo. Las condiciones de borde tomadas son las correspondientes a la viga empotrada en un extremo, 
donde la tangente y el desplazamiento en ese punto son nulos. 




Vista seccion media 
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d 2 y = m(r) 
dR 2 E-J 

P = 4 • n • p • (k + ^/l + tan'p)- sin(a + p) 

O(R min )=0 ^ I 3 (R min ) = -C 3 

I 3 (R) = ^R 6 + ^R 4 + ^R 5 + ^R 3 + ^R 2 + C 2 R 
120 24 60 6 2 

y(R) = Jo(R)dR = ^ y Jl 3 (R)dR = ^ y [l 4 (R) + C 4 ] 
y(R mi „) = o ^ I 4 (R min ) = -C 4 

I 4 (R) = ^R 7 + ^R 5 + ^R 6 + ^R 4 + ^R 3 + ^R 2 + C 3 R 
v y 840 120 360 24 6 2 

La deflexion maxima se dara en el extremo de la pala libre (R max ) Para el momenta de inercia de la 
seccion media y el modulo de Young del acero 1020 (200 GPa) la deflexion sera: 

E = 200GPa = 2-10 n Pa 
J m edio= 0.00180018 

y(R max ) = -^-[l 4 (R max ) + C 4 ] = 3.539m 
E -J 

Que teniendo en cuenta que es la deflexion tomada para el peor de los casos puede considerarse la 

deflexion maxima que tendra la pala en condiciones normales de funcionamiento. 

Analisis estatico de la columna 

Tensiones de flexion debido al empuje de las palas 

Se tomara una columna de seccion anular, teniendo en cuenta que para ascender a la cima de la torre 
por dentro de la columna debe correr una escalera o ascensor. Su altura debe ser al menos del radio de la 
pala mas una distancia de seguridad. La mayoria de las turbinas consultadas tienen la columna del mismo 
tamano que el diametro de barrido, como la turbina del presente trabajo es de dimensiones mucho mayores, 
se toma en primera aproximacion el alto de la columna como una vez y media del radio maximo de las 
palas. La altura de la columna debera ser tal que el viento incidente en el tubo que debe desplazar la turbina 
sea mas o menos constante en su area. 

En este caso se tomara como fuerza predominante la de las palas debido a que es mucho mayor a las 
demas. Para eso se debe integrar la fuerza distribuida sobre las mismas. La integral definida de la fuerza 
entre R min y R ma x coincide con la funcion del esfuerzo cortante evaluada entre esos mismos puntos de 
manera que la fuerza total en las palas sera tres veces ese valor. 
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4 • 7i • p • (bR + c)- (k + Vl+tan'pj- V a 2 • sin(a + p) 
= 3 • 4 • tc • p • (k + ^/l+tan'p)- sin(a + p) • J (bR + c)((coR) 2 + V v 2 )dR 



F total =12-TC-p-(K + Vl+tan 2 p)-sin(a + p)-[l 1 (R max )-I 1 (R min )] 
I x (r) = j"(bR 3 co 2 + bV v 2 R + cR V + cV v 2 )• dR = 

I 1 (R) = ^R 4 + ^R 2 + ^R 3 + cV v 2 R 
V ' 4 2 3 

F total = 3.4488 -10 5 = 344.88 KN 

En este caso la fuerza se tomara como concentrada en un extreme De forma que puede calcularse 
facilmente el momenta flector en la base (la seccion mas solicitada) y calcular la tension maxima. 
M 



J = ^( D 4 -d 4 )=^D 



64 



14 \ 



1- 



D z 



— D 4 (l-5 4 ) 
64 v ' 



64 

D 
_ D 

J max ,-v 

M = F total -L 

a = 32 ' F -'' L 
max rcD 3 (l-5 4 ) 

Tensiones axiales debido al peso de la torre 

Ademas de las tensiones de flexion deben tomarse en cuenta las tensiones axiales debido al peso de 
la maquinaria y las palas. El peso de cada pala es sacado del programa de CAD. Dado que el peso 
predominante es el de las palas, se toma sin mayor error el triple del peso de una pala (aproximadamente 
300 toneladas). 

Las tensiones axiales seran entonces: 
P 



a = 



A=*( D 2 -d 2 )=*D ; 



8 = 



4 

d_ 

D 



|2 A 



D^ 



= *D 2 (l-S 2 ) 
4 v ; 



P = 3-102451.36 Kg • 9.81!% = 3.015M0 6 N = 3.0151 MN 



a = 



4P 

ttD 2 (1-5 2 ) 



Dimensionamiento de la seccion de la columna 

Las tensiones totales seran la suma de ambas tensiones. 
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° flexion + ° axial 



4-P 32-F total -L + 4-P-D-(l+S 2 ) 

7rD 3 (l-8 4 )^7rD 2 (l-8 2 ) 



32-F total -L 



ttD 3 (1-5 4 ) 

Tomando la relacion y 8' y fijado el limite de tension conociendo el limite de fluencia del material 
afectado por un coeficiente de seguridad es posible determinar el diametro exterior. 

En nuestro caso se adoptara un diametro interior a fin de asegurar el problema de la accesibilidad a la 
torre desde el interior y con ese diametro se fijara el exterior. Adoptaremos un diametro interior de 3 metros 
y un diametro exterior de 3.3 metros y se verificara la tension resultante comparandola con la tension limite 
elastica del acero 1020. 

D = 3.3m 



d = 3m 
3m 



= 0.9091 



3.3m 
a„ = 351.571 MPa 



CT 



32 • F total • L 



4-P 



ttD 3 (1-8 4 ) 7iD 2 (l-8 2 ) 



2.4237 -10 7 = 24.237MPa < a f 



Analisis dinamico de la pala 

Frecuencias naturales de vibracion en flexion de la pala 

Es conocida la ecuacion que debe resolverse para obtener las frecuencias de vibracion natural de la 
viga en flexion. En el caso de las palas se las modeliza como una viga empotrada en un extremo, 
considerando un momenta de inercia de la seccion constante e igual al momenta de inercia de la seccion 
media. 

_i 

cos(u)-cosh(u) + l= => cos(u) = 1-^- 

cosh(u) 




Considerando 'a' como constante, se puede resolver la ecuacion numericamente para V para despues 
hallar los valores de V. Para tener alguna idea de los valores de V para los que se cumple la ecuacion, se 
graficaran ambos miembros en un mismo grafico. 





cos(x) 
-l/cosh(x) 


1 

0.5 



0.5 

-1 

1 5 




___^__;__;£_^V- H L - _\ 1 __/_ 1 _\_ + __ + _A_\ / - A ! 1- _ _/ _ _\ 

\ ' '/ \ / \ / \ / \ / V ' /' 


■ \ i /: \ / rN K7 \ / i\ i / i - 


! ! 1 1 1 - 


ii iii 
ii i i i 



-20 -18 -16 -14 -12 -10 -8 -6 -4 -2 2 4 6 8 10 12 14 16 18 20 



Grafico cos(x)--l/cosh(x) 
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De lo que podemos ver que las curvas se cruzan entre 1 y 2, entre 4 y 5 y luego tienden a los ceros 
del coseno trigonometrico. Resolviendo las primeras frecuencias de forma numerica se obtuvo: 



Uj =1.8751 



u 2 =4.6941 

u 3 = 7.8548 

Lo que resulta en las siguientes frecuencias respectivas: 




CO : 



E-J 



media 



J = J 



medio 



0.276803 m' 
0.00180018 m 4 



p = 7900 
l = 48m 



Kg/ 



m 



CO, 



0.6192 V 



'u> 2 



W 



=> a = 



/E-J 
^A 



= 405.7645 



co, = 3.8806 



co, =10.8657 y 

3 A 



Siendo la velocidad de rotacion de la turbina: 



co r 



1.5708 V 



La turbina disenada funciona con una frecuencia de giro V situada entre la primera y la segunda 
frecuencia de resonancia. Debido a que la primera frecuencia de resonancia es muy baja, es posible evitarla, 
al parar la turbina, frenandola lo suficientemente rapido como para que el paso por esa frecuencia no la 
Neve a falla. 



Analisis dinamico de la columna 

En la columna es necesario efectuar dos analisis de frecuencias naturales, el de vibraciones axiales y el 
de vibraciones de flexion, dado que ella debe soportar el peso de toda la estructura asi como el empuje de 
la turbina que la llevara a flexionarse. 

Frecuencias naturales de vibracion axiales 

En este caso es necesario tomar las vibraciones axiales de una barra con una masa concentrada en el 
extremo, cuya ecuacion es: 
B 



a 



m 



tan (a) 



a = VX -1 

Donde 'B' es la masa de la barra y 'm' es la masa concentrada, que tomaremos como el triple del peso 
de una pala (cerca de 300 toneladas) considerando el resto de las masas despreciables. 
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p = 7900 K g/ 3 
/m 

1 = 48m -1.5 = 72 m 

A = -(D 2 -d 2 ) = 1.4844 m 2 
4 V ; 

B = p-1-A = 846.60192 ton 

m = 3 • 102.45136 ton = 307.35408 ton 

B 



m 



a • tan(a) 



a = v X • 1 



CO 



Xc- 



a 



1 VP 
Nuevamente se graficaran ambos miembros de la ecuacion para determinar los valores de 'a' y luego 



los de las frecuencias naturales. 



x tan(x) 
B/m 



60 



40 



20 



-20 



-40 



-60 



I _ L _ _ 1 



.1 I I I 



-20 -18 -16 -14 -12 -10 -8 -6 -4 



2 2 4 6 

x 

Grafico x*tan(x)-B/m 



8 10 12 14 16 18 20 



Donde se pueden obtener numericamente los valores de V. 
a l =1.1693 

a 2 =1.5708 
a 3 = 3.7723 

Y sus correspondientes valores de frecuencia. 
a IE 

Hp 



CO 



co, = 81.7137 V 

1 /s 

co, =109.7715 y 

2 /s 

co, =263.6181 y 

3 /s 

co rot =1.5708l{ 

Contrastados con la frecuencia de rotacion del rotor se puede asegurar de estar suficientemente 
alejados de cualquier situacion de resonancia. 
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Vibraciones de flexion 

Nuevamente la ecuacion a resolver es 
cos(u)-cosh(u) + l= => cos(u) : 



cosh(u) 



•1 



E-J 
f^A 



CO : 



E-J 

p^A 



^u^ 2 



W 



Dado que el parametro V de la ecuacion de frecuencias naturales de vibracion en flexion es 
adimensional, es posible utilizar los mismos valores para determinar las frecuencias de resonancia sobre la 
columna en flexion. Cambian el parametro Vy la longitud de la columna. 

E = 200GPa 



A = -(D 2 -d 2 ) = 1.4844 



m' 



j = — (D 4 -d 4 ) =1.845298 m' 

64 v ' 



p = 7900 
l = 72m 



K §/ 



m 



E-J 
>A 



5609.95316 



Uj =1.8751 
u 2 =4.6941 
u 3 = 7.8548 
co, =3.80490 



co 2 = 23.84509 

co, = 66.76740 V 

3 /s 

Y la frecuencia del rotor 

co rot = 1.5708 y s 

Que se encuentra muy por debajo de la frecuencia de resonancia de la columna. 
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26/12/08 02:49 C:\Documents and Settings\Joaquin\Escritorio . . AGraf icoperf il ,m 1 of 1 

function Graficoperfil(K, alfad, beta, Rmin, Rmax, amin, amax, wd, Vd) 

% Funcion que dibuja varios perfiles para generar una vista preeliminar en 
% 3d de la pala 

b= ( amin - amax ) / ( Rmax - Rmin ) ; 
c=amax-b*Rmin; 

for R=Rmin:Rmax 

a=b*R+c; 

lambdad=atan ( wd*R/Vd ) ; 

perf il(K, a, alfad, beta, lambdad, R) ; 
end 
perf il3d(Rmin, Rmax) ; 



26/12/08 02:48 C:\Documents and Settings\Joaquin\Escritorio\Perfil\perf il.m 1 of 1 

function perfil(K, a, alfad, beta, lambdad, R) 

% Dibuja la curva del contorno del perfil 

% aerodinamico segun la transformacion de Kutta- Jukovsky 

% dadas las variables: 

% K = Relacion de d/a para calcular las dimensiones del perfil 

% a = Puntos de corte de la circunf erencia del esqueleto del perfil 

% alfad = angulo de ataque (solo usado para el giro del perfil) 

% beta = angulo de arqueo del perfil 

% lambdad = angulo de rotacion del perfil (usado para el giro del perfil y generar alas * 

en 3d) 

% R = posicion r (usado para generar alas en 3d) 



%% Generacion de puntos para el perfil 

d=a*K; 

r=a/cos(beta) ; 

m=r*sin(beta) ; 

R0=r+d ; 

t=[(-pi+beta) :2*pi/36: (pi+beta)] ' ; 

xa=r*cos(t ) ; ya=r*sin(t) ; 

za=xa+i*(ya+m) ; 

wa=( (za+a A 2 . /za) -i*2*m) *exp(i*( -pi/2+lambdad+alf ad) ) ; 

xp=R0*cos(t) ; yp=R0*sin(t) ; 

zper=xp+i*yp+d *cos ( beta )+i*(m+d* sin (beta) ) ; 

wp=( (zper+a A 2./zper) -i*2*m)*exp(i*( -pi/2+lambdad+alfad) ) ; 

%% Salida de datos en texto para Solidworks 

E=1000*[real(wp) imag(wp)]; 

E( :,3)=R*1000; 

filename=[ ' perfilR' num2str(R) '.txt']; 

dlmwrite( filename, E, ' delimiter ' , ' \t ' , ' precision ' ,4, ' newline ' , ' pc ' ) ; 

for k=l:length(t)/2 

E=1000*[real(wa) imag(wa)]; 

E( :,3)=R*1000; 

filename=[ 'Almat ' num2str(k) '.txt']; 

dlmwrite( filename, E(k, : ), ' -append ' , 'delimiter ' , ' \t ' , ' precision ' ,4, ' newline' , ' pc ' ) 
end 

%% Grafico de perfil 

Perfil=figure( ' name' ,[' Perfil R = ' num2str(R)] ) ; 

ejes=axes( . . . 

'Visible', 'on', . . . 

' Box' , 'on ' , . . . 

' DataAspectRatioMode ' , ' manual ' , . . . 

' DataAspectRatio' , [1 1 1] , . . . 

' Parent ' , Perfil) ; 
xlabel( 'x [mm] ' ) ; 
ylabel( 'y [mm] ' ) ; 
hold on 

plot ( real (wp)*1000,imag(wp)*1000, 'Parent' ,ejes) ; 
plot(real(wa)*1000,imag(wa)*1000, ' -r ' , 'Parent' ,ejes); 
hold off 
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Mass properties of Pala3 ( Part Configuration - Default ) 

Output coordinate System: - default -- 

Density = 7900 kilograms per cubic meter 

Mass = 102451.36 kilograms 

Volume = 12.968527 cubic meters 

Surface area = 139.83828 meters A 2 

Center of mass: ( meters ) 
X = 0.29393554 
Y = -0.14322289 
Z = 19.886729 

Principal axes of inertia and principal moments of inertia: ( kilograms * square meters ) 

Taken at the center of mass. 
Ix = (-0.0035066039, 0.0077259513, 0.99996401) Px = 22449.766 
Iy = (0.89537773, -0.44525882, 0.0065800199) Py = 12281506 
Iz = (0.44529363, 0.89536858, -0.0053562984) Pz = 12299873 

Moments of inertia: ( kilograms * square meters ) 

Taken at the center of mass and aligned with the output coordinate system. 

Lxx = 12284997 Lxy = -7655.1791 Lxz = -42942.299 

Lyx = -7655.1791 Lyy = 12295499 Lyz = 94797.548 

Lzx = -42942.299 Lzy = 94797.548 Lzz = 23332.781 

Moments of inertia: ( kilograms * square meters ) 
Taken at the output coordinate system. 

Ixx = 52804767 lxy = -11968.207 lxz = 555928.57 

lyx = -11968.207 lyy = 52822019 lyz = -197007.98 

lzx = 555928.57 lzy = -197007.98 lzz = 34285.948 
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Pala, Seccion R=5m 






Measurements are based on sectioned model 
Section properties of the selected faces of Pala3 

Area = 0.62264829 meters A 2 

Centroid relative to output coordinate system origin: ( meters ) 
X = 0.2605922 
Y = -0.35013616 
Z = 5.011 

Moments of inertia of the area, at the centroid: ( meters A 4 ) 
Lxx = 0.13211298 Lxy = 0.09659872 Lxz = 
Lyx = 0.09659872 Lyy = 0.084970669 Lyz = 
Lzx = Lzy = Lzz = 0.21708365 



Polar moment of inertia of the area, at the centroid = 0.21708365 meters A 4 

Angle between principal axes and part axes = 128.14358 degrees 

Principal moments of inertia of the area, at the centroid: ( meters A 4 ) 
Ix = 0.0091088732 
Iy = 0.20797478 



Moments of inertia of the area, at the output coordinate system: ( meters A 4 ) 
LXX = 15.843221 LXY = -0.15341086 LXZ = 0.81307126 
LYX = -0.15341086 LYY = 15.762027 LYZ = -1.0924565 
LZX = 0.81307126 LZY = -1.0924565 LZZ = 0.33570041 




R 255,973 mm 




T 



Pala, Seccion media 



Measurements are based on sectioned model 
Section properties of the selected faces of Pala3 

Area = 0.27680347 meters A 2 

Centroid relative to output coordinate system origin: ( meters ) 
X = 0.28286278 
Y = -0.067887056 
Z = 26.5 

Moments of inertia of the area, at the centroid: ( meters A 4 ) 
Lxx = 0.0035037569 Lxy = 0.0080034882 Lxz = 
Lyx = 0.0080034882 Lyy = 0.039401072 Lyz = 
Lzx = Lzy = Lzz = 0.042904829 

Polar moment of inertia of the area, at the centroid = 0.042904829 meters A 4 

Angle between principal axes and part axes = -12.016288 degrees 

Principal moments of inertia of the area, at the centroid: ( meters A 4 ) 
Ix = 0.0018001848 
Iy = 0.041104644 



Moments of inertia of the area, at the output coordinate system: ( meters A 4 ' 
LXX = 194.39001 LXY = -0.013318868 LXZ = 2.0748811 
LYX = -0.013318868 LYY = 194.44678 LYZ = -0.49797137 
LZX = 2.0748811 LZY = -0.49797137 LZZ = 0.066327941 
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